3.2.2 \(\int \frac {x^3 \sin (c+d x)}{(a+b x^3)^2} \, dx\) [102]

3.2.2.1 Optimal result
3.2.2.2 Mathematica [C] (verified)
3.2.2.3 Rubi [A] (verified)
3.2.2.4 Maple [C] (verified)
3.2.2.5 Fricas [C] (verification not implemented)
3.2.2.6 Sympy [F(-1)]
3.2.2.7 Maxima [F]
3.2.2.8 Giac [F]
3.2.2.9 Mupad [F(-1)]

3.2.2.1 Optimal result

Integrand size = 19, antiderivative size = 714 \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^3\right )^2} \, dx=-\frac {(-1)^{2/3} d \cos \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 \sqrt [3]{a} b^{5/3}}-\frac {d \cos \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\sqrt [3]{-1} d \cos \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\operatorname {CosIntegral}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right ) \sin \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 a^{2/3} b^{4/3}}-\frac {\sqrt [3]{-1} \operatorname {CosIntegral}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right ) \sin \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 a^{2/3} b^{4/3}}+\frac {(-1)^{2/3} \operatorname {CosIntegral}\left (\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right ) \sin \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{9 a^{2/3} b^{4/3}}-\frac {x \sin (c+d x)}{3 b \left (a+b x^3\right )}+\frac {\sqrt [3]{-1} \cos \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 a^{2/3} b^{4/3}}-\frac {(-1)^{2/3} d \sin \left (c+\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {\cos \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 a^{2/3} b^{4/3}}+\frac {d \sin \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 \sqrt [3]{a} b^{5/3}}+\frac {(-1)^{2/3} \cos \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 a^{2/3} b^{4/3}}-\frac {\sqrt [3]{-1} d \sin \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}+d x\right )}{9 \sqrt [3]{a} b^{5/3}} \]

output
-1/9*d*Ci(a^(1/3)*d/b^(1/3)+d*x)*cos(c-a^(1/3)*d/b^(1/3))/a^(1/3)/b^(5/3)- 
1/9*(-1)^(2/3)*d*Ci((-1)^(1/3)*a^(1/3)*d/b^(1/3)-d*x)*cos(c+(-1)^(1/3)*a^( 
1/3)*d/b^(1/3))/a^(1/3)/b^(5/3)+1/9*(-1)^(1/3)*d*Ci((-1)^(2/3)*a^(1/3)*d/b 
^(1/3)+d*x)*cos(c-(-1)^(2/3)*a^(1/3)*d/b^(1/3))/a^(1/3)/b^(5/3)-1/9*(-1)^( 
1/3)*cos(c+(-1)^(1/3)*a^(1/3)*d/b^(1/3))*Si(-(-1)^(1/3)*a^(1/3)*d/b^(1/3)+ 
d*x)/a^(2/3)/b^(4/3)+1/9*cos(c-a^(1/3)*d/b^(1/3))*Si(a^(1/3)*d/b^(1/3)+d*x 
)/a^(2/3)/b^(4/3)+1/9*(-1)^(2/3)*cos(c-(-1)^(2/3)*a^(1/3)*d/b^(1/3))*Si((- 
1)^(2/3)*a^(1/3)*d/b^(1/3)+d*x)/a^(2/3)/b^(4/3)+1/9*Ci(a^(1/3)*d/b^(1/3)+d 
*x)*sin(c-a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)+1/9*d*Si(a^(1/3)*d/b^(1/3)+d* 
x)*sin(c-a^(1/3)*d/b^(1/3))/a^(1/3)/b^(5/3)-1/9*(-1)^(1/3)*Ci((-1)^(1/3)*a 
^(1/3)*d/b^(1/3)-d*x)*sin(c+(-1)^(1/3)*a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)+ 
1/9*(-1)^(2/3)*d*Si(-(-1)^(1/3)*a^(1/3)*d/b^(1/3)+d*x)*sin(c+(-1)^(1/3)*a^ 
(1/3)*d/b^(1/3))/a^(1/3)/b^(5/3)+1/9*(-1)^(2/3)*Ci((-1)^(2/3)*a^(1/3)*d/b^ 
(1/3)+d*x)*sin(c-(-1)^(2/3)*a^(1/3)*d/b^(1/3))/a^(2/3)/b^(4/3)-1/9*(-1)^(1 
/3)*d*Si((-1)^(2/3)*a^(1/3)*d/b^(1/3)+d*x)*sin(c-(-1)^(2/3)*a^(1/3)*d/b^(1 
/3))/a^(1/3)/b^(5/3)-1/3*x*sin(d*x+c)/b/(b*x^3+a)
 
3.2.2.2 Mathematica [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 4 in optimal.

Time = 0.17 (sec) , antiderivative size = 383, normalized size of antiderivative = 0.54 \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^3\right )^2} \, dx=\frac {\text {RootSum}\left [a+b \text {$\#$1}^3\&,\frac {i \cos (c+d \text {$\#$1}) \operatorname {CosIntegral}(d (x-\text {$\#$1}))+\operatorname {CosIntegral}(d (x-\text {$\#$1})) \sin (c+d \text {$\#$1})+\cos (c+d \text {$\#$1}) \text {Si}(d (x-\text {$\#$1}))-i \sin (c+d \text {$\#$1}) \text {Si}(d (x-\text {$\#$1}))+d \cos (c+d \text {$\#$1}) \operatorname {CosIntegral}(d (x-\text {$\#$1})) \text {$\#$1}-i d \operatorname {CosIntegral}(d (x-\text {$\#$1})) \sin (c+d \text {$\#$1}) \text {$\#$1}-i d \cos (c+d \text {$\#$1}) \text {Si}(d (x-\text {$\#$1})) \text {$\#$1}-d \sin (c+d \text {$\#$1}) \text {Si}(d (x-\text {$\#$1})) \text {$\#$1}}{\text {$\#$1}^2}\&\right ]+\text {RootSum}\left [a+b \text {$\#$1}^3\&,\frac {-i \cos (c+d \text {$\#$1}) \operatorname {CosIntegral}(d (x-\text {$\#$1}))+\operatorname {CosIntegral}(d (x-\text {$\#$1})) \sin (c+d \text {$\#$1})+\cos (c+d \text {$\#$1}) \text {Si}(d (x-\text {$\#$1}))+i \sin (c+d \text {$\#$1}) \text {Si}(d (x-\text {$\#$1}))+d \cos (c+d \text {$\#$1}) \operatorname {CosIntegral}(d (x-\text {$\#$1})) \text {$\#$1}+i d \operatorname {CosIntegral}(d (x-\text {$\#$1})) \sin (c+d \text {$\#$1}) \text {$\#$1}+i d \cos (c+d \text {$\#$1}) \text {Si}(d (x-\text {$\#$1})) \text {$\#$1}-d \sin (c+d \text {$\#$1}) \text {Si}(d (x-\text {$\#$1})) \text {$\#$1}}{\text {$\#$1}^2}\&\right ]-\frac {6 b x \sin (c+d x)}{a+b x^3}}{18 b^2} \]

input
Integrate[(x^3*Sin[c + d*x])/(a + b*x^3)^2,x]
 
output
(RootSum[a + b*#1^3 & , (I*Cos[c + d*#1]*CosIntegral[d*(x - #1)] + CosInte 
gral[d*(x - #1)]*Sin[c + d*#1] + Cos[c + d*#1]*SinIntegral[d*(x - #1)] - I 
*Sin[c + d*#1]*SinIntegral[d*(x - #1)] + d*Cos[c + d*#1]*CosIntegral[d*(x 
- #1)]*#1 - I*d*CosIntegral[d*(x - #1)]*Sin[c + d*#1]*#1 - I*d*Cos[c + d*# 
1]*SinIntegral[d*(x - #1)]*#1 - d*Sin[c + d*#1]*SinIntegral[d*(x - #1)]*#1 
)/#1^2 & ] + RootSum[a + b*#1^3 & , ((-I)*Cos[c + d*#1]*CosIntegral[d*(x - 
 #1)] + CosIntegral[d*(x - #1)]*Sin[c + d*#1] + Cos[c + d*#1]*SinIntegral[ 
d*(x - #1)] + I*Sin[c + d*#1]*SinIntegral[d*(x - #1)] + d*Cos[c + d*#1]*Co 
sIntegral[d*(x - #1)]*#1 + I*d*CosIntegral[d*(x - #1)]*Sin[c + d*#1]*#1 + 
I*d*Cos[c + d*#1]*SinIntegral[d*(x - #1)]*#1 - d*Sin[c + d*#1]*SinIntegral 
[d*(x - #1)]*#1)/#1^2 & ] - (6*b*x*Sin[c + d*x])/(a + b*x^3))/(18*b^2)
 
3.2.2.3 Rubi [A] (verified)

Time = 1.51 (sec) , antiderivative size = 725, normalized size of antiderivative = 1.02, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.263, Rules used = {3824, 3814, 2009, 3827, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {x^3 \sin (c+d x)}{\left (a+b x^3\right )^2} \, dx\)

\(\Big \downarrow \) 3824

\(\displaystyle \frac {\int \frac {\sin (c+d x)}{b x^3+a}dx}{3 b}+\frac {d \int \frac {x \cos (c+d x)}{b x^3+a}dx}{3 b}-\frac {x \sin (c+d x)}{3 b \left (a+b x^3\right )}\)

\(\Big \downarrow \) 3814

\(\displaystyle \frac {\int \left (-\frac {\sin (c+d x)}{3 a^{2/3} \left (-\sqrt [3]{b} x-\sqrt [3]{a}\right )}-\frac {\sin (c+d x)}{3 a^{2/3} \left (\sqrt [3]{-1} \sqrt [3]{b} x-\sqrt [3]{a}\right )}-\frac {\sin (c+d x)}{3 a^{2/3} \left (-(-1)^{2/3} \sqrt [3]{b} x-\sqrt [3]{a}\right )}\right )dx}{3 b}+\frac {d \int \frac {x \cos (c+d x)}{b x^3+a}dx}{3 b}-\frac {x \sin (c+d x)}{3 b \left (a+b x^3\right )}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {d \int \frac {x \cos (c+d x)}{b x^3+a}dx}{3 b}+\frac {\frac {\sin \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}-\frac {\sqrt [3]{-1} \sin \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \operatorname {CosIntegral}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {(-1)^{2/3} \sin \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {\sqrt [3]{-1} \cos \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \text {Si}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {\cos \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {(-1)^{2/3} \cos \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}}{3 b}-\frac {x \sin (c+d x)}{3 b \left (a+b x^3\right )}\)

\(\Big \downarrow \) 3827

\(\displaystyle \frac {d \int \left (-\frac {\cos (c+d x)}{3 \sqrt [3]{a} \sqrt [3]{b} \left (\sqrt [3]{b} x+\sqrt [3]{a}\right )}-\frac {(-1)^{2/3} \cos (c+d x)}{3 \sqrt [3]{a} \sqrt [3]{b} \left (\sqrt [3]{a}-\sqrt [3]{-1} \sqrt [3]{b} x\right )}+\frac {\sqrt [3]{-1} \cos (c+d x)}{3 \sqrt [3]{a} \sqrt [3]{b} \left ((-1)^{2/3} \sqrt [3]{b} x+\sqrt [3]{a}\right )}\right )dx}{3 b}+\frac {\frac {\sin \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}-\frac {\sqrt [3]{-1} \sin \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \operatorname {CosIntegral}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {(-1)^{2/3} \sin \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {\sqrt [3]{-1} \cos \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \text {Si}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {\cos \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {(-1)^{2/3} \cos \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}}{3 b}-\frac {x \sin (c+d x)}{3 b \left (a+b x^3\right )}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {\frac {\sin \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}-\frac {\sqrt [3]{-1} \sin \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \operatorname {CosIntegral}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {(-1)^{2/3} \sin \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {\sqrt [3]{-1} \cos \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \text {Si}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {\cos \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}+\frac {(-1)^{2/3} \cos \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 a^{2/3} \sqrt [3]{b}}}{3 b}+\frac {d \left (-\frac {(-1)^{2/3} \cos \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \operatorname {CosIntegral}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{3 \sqrt [3]{a} b^{2/3}}-\frac {\cos \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 \sqrt [3]{a} b^{2/3}}+\frac {\sqrt [3]{-1} \cos \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \operatorname {CosIntegral}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 \sqrt [3]{a} b^{2/3}}-\frac {(-1)^{2/3} \sin \left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}+c\right ) \text {Si}\left (\frac {\sqrt [3]{-1} \sqrt [3]{a} d}{\sqrt [3]{b}}-d x\right )}{3 \sqrt [3]{a} b^{2/3}}+\frac {\sin \left (c-\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (x d+\frac {\sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 \sqrt [3]{a} b^{2/3}}-\frac {\sqrt [3]{-1} \sin \left (c-\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right ) \text {Si}\left (x d+\frac {(-1)^{2/3} \sqrt [3]{a} d}{\sqrt [3]{b}}\right )}{3 \sqrt [3]{a} b^{2/3}}\right )}{3 b}-\frac {x \sin (c+d x)}{3 b \left (a+b x^3\right )}\)

input
Int[(x^3*Sin[c + d*x])/(a + b*x^3)^2,x]
 
output
-1/3*(x*Sin[c + d*x])/(b*(a + b*x^3)) + ((CosIntegral[(a^(1/3)*d)/b^(1/3) 
+ d*x]*Sin[c - (a^(1/3)*d)/b^(1/3)])/(3*a^(2/3)*b^(1/3)) - ((-1)^(1/3)*Cos 
Integral[((-1)^(1/3)*a^(1/3)*d)/b^(1/3) - d*x]*Sin[c + ((-1)^(1/3)*a^(1/3) 
*d)/b^(1/3)])/(3*a^(2/3)*b^(1/3)) + ((-1)^(2/3)*CosIntegral[((-1)^(2/3)*a^ 
(1/3)*d)/b^(1/3) + d*x]*Sin[c - ((-1)^(2/3)*a^(1/3)*d)/b^(1/3)])/(3*a^(2/3 
)*b^(1/3)) + ((-1)^(1/3)*Cos[c + ((-1)^(1/3)*a^(1/3)*d)/b^(1/3)]*SinIntegr 
al[((-1)^(1/3)*a^(1/3)*d)/b^(1/3) - d*x])/(3*a^(2/3)*b^(1/3)) + (Cos[c - ( 
a^(1/3)*d)/b^(1/3)]*SinIntegral[(a^(1/3)*d)/b^(1/3) + d*x])/(3*a^(2/3)*b^( 
1/3)) + ((-1)^(2/3)*Cos[c - ((-1)^(2/3)*a^(1/3)*d)/b^(1/3)]*SinIntegral[(( 
-1)^(2/3)*a^(1/3)*d)/b^(1/3) + d*x])/(3*a^(2/3)*b^(1/3)))/(3*b) + (d*(-1/3 
*((-1)^(2/3)*Cos[c + ((-1)^(1/3)*a^(1/3)*d)/b^(1/3)]*CosIntegral[((-1)^(1/ 
3)*a^(1/3)*d)/b^(1/3) - d*x])/(a^(1/3)*b^(2/3)) - (Cos[c - (a^(1/3)*d)/b^( 
1/3)]*CosIntegral[(a^(1/3)*d)/b^(1/3) + d*x])/(3*a^(1/3)*b^(2/3)) + ((-1)^ 
(1/3)*Cos[c - ((-1)^(2/3)*a^(1/3)*d)/b^(1/3)]*CosIntegral[((-1)^(2/3)*a^(1 
/3)*d)/b^(1/3) + d*x])/(3*a^(1/3)*b^(2/3)) - ((-1)^(2/3)*Sin[c + ((-1)^(1/ 
3)*a^(1/3)*d)/b^(1/3)]*SinIntegral[((-1)^(1/3)*a^(1/3)*d)/b^(1/3) - d*x])/ 
(3*a^(1/3)*b^(2/3)) + (Sin[c - (a^(1/3)*d)/b^(1/3)]*SinIntegral[(a^(1/3)*d 
)/b^(1/3) + d*x])/(3*a^(1/3)*b^(2/3)) - ((-1)^(1/3)*Sin[c - ((-1)^(2/3)*a^ 
(1/3)*d)/b^(1/3)]*SinIntegral[((-1)^(2/3)*a^(1/3)*d)/b^(1/3) + d*x])/(3*a^ 
(1/3)*b^(2/3))))/(3*b)
 

3.2.2.3.1 Defintions of rubi rules used

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 3814
Int[((a_) + (b_.)*(x_)^(n_))^(p_)*Sin[(c_.) + (d_.)*(x_)], x_Symbol] :> Int 
[ExpandIntegrand[Sin[c + d*x], (a + b*x^n)^p, x], x] /; FreeQ[{a, b, c, d}, 
 x] && ILtQ[p, 0] && IGtQ[n, 0] && (EqQ[n, 2] || EqQ[p, -1])
 

rule 3824
Int[(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_)*Sin[(c_.) + (d_.)*(x_)], x_Sym 
bol] :> Simp[x^(m - n + 1)*(a + b*x^n)^(p + 1)*(Sin[c + d*x]/(b*n*(p + 1))) 
, x] + (-Simp[(m - n + 1)/(b*n*(p + 1))   Int[x^(m - n)*(a + b*x^n)^(p + 1) 
*Sin[c + d*x], x], x] - Simp[d/(b*n*(p + 1))   Int[x^(m - n + 1)*(a + b*x^n 
)^(p + 1)*Cos[c + d*x], x], x]) /; FreeQ[{a, b, c, d, m}, x] && ILtQ[p, -1] 
 && IGtQ[n, 0] && (GtQ[m - n + 1, 0] || GtQ[n, 2]) && RationalQ[m]
 

rule 3827
Int[Cos[(c_.) + (d_.)*(x_)]*(x_)^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_), x_Sym 
bol] :> Int[ExpandIntegrand[Cos[c + d*x], x^m*(a + b*x^n)^p, x], x] /; Free 
Q[{a, b, c, d, m}, x] && ILtQ[p, 0] && IGtQ[n, 0] && (EqQ[n, 2] || EqQ[p, - 
1]) && IntegerQ[m]
 
3.2.2.4 Maple [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.93 (sec) , antiderivative size = 1184, normalized size of antiderivative = 1.66

method result size
derivativedivides \(\text {Expression too large to display}\) \(1184\)
default \(\text {Expression too large to display}\) \(1184\)
risch \(\text {Expression too large to display}\) \(1379\)

input
int(x^3*sin(d*x+c)/(b*x^3+a)^2,x,method=_RETURNVERBOSE)
 
output
1/d^4*(-d^6*c^3*(sin(d*x+c)*(1/3/a/d^3*(d*x+c)-1/3*c/a/d^3)/(a*d^3-c^3*b+3 
*b*c^2*(d*x+c)-3*b*c*(d*x+c)^2+b*(d*x+c)^3)+2/9/a/d^3/b*sum(1/(_R1^2-2*_R1 
*c+c^2)*(-Si(-d*x+_R1-c)*cos(_R1)+Ci(d*x-_R1+c)*sin(_R1)),_R1=RootOf(_Z^3* 
b-3*_Z^2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))+1/9/a/d^3/b*sum(1/(-_RR1+c)*(Si(-d*x 
+_RR1-c)*sin(_RR1)+Ci(d*x-_RR1+c)*cos(_RR1)),_RR1=RootOf(_Z^3*b-3*_Z^2*b*c 
+3*_Z*b*c^2+a*d^3-b*c^3)))+sin(d*x+c)*(c^2*d^3/a*(d*x+c)^2-c^3*d^3/a*(d*x+ 
c))/(a*d^3-c^3*b+3*b*c^2*(d*x+c)-3*b*c*(d*x+c)^2+b*(d*x+c)^3)+1/3*c^2*d^3/ 
a/b*sum((c+_R1)/(_R1^2-2*_R1*c+c^2)*(-Si(-d*x+_R1-c)*cos(_R1)+Ci(d*x-_R1+c 
)*sin(_R1)),_R1=RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))+1/3*c^2* 
d^3/a/b*sum(_RR1/(-_RR1+c)*(Si(-d*x+_RR1-c)*sin(_RR1)+Ci(d*x-_RR1+c)*cos(_ 
RR1)),_RR1=RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))+sin(d*x+c)*(- 
2*c^2*d^3/a*(d*x+c)^2+3*c^3*d^3/a*(d*x+c)+c*d^3*(a*d^3-b*c^3)/a/b)/(a*d^3- 
c^3*b+3*b*c^2*(d*x+c)-3*b*c*(d*x+c)^2+b*(d*x+c)^3)-2/3*c^2*d^3/a/b*sum(_R1 
/(_R1^2-2*_R1*c+c^2)*(-Si(-d*x+_R1-c)*cos(_R1)+Ci(d*x-_R1+c)*sin(_R1)),_R1 
=RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z*b*c^2+a*d^3-b*c^3))-1/3*c*d^3/a/b^2*sum((-2 
*_RR1^2*b*c+3*_RR1*b*c^2+a*d^3-b*c^3)/(_RR1^2-2*_RR1*c+c^2)*(Si(-d*x+_RR1- 
c)*sin(_RR1)+Ci(d*x-_RR1+c)*cos(_RR1)),_RR1=RootOf(_Z^3*b-3*_Z^2*b*c+3*_Z* 
b*c^2+a*d^3-b*c^3))+sin(d*x+c)*(c^2*d^3/a*(d*x+c)^2-1/3*d^3*(a*d^3+5*b*c^3 
)/a/b*(d*x+c)-2/3*c*d^3*(a*d^3-b*c^3)/a/b)/(a*d^3-c^3*b+3*b*c^2*(d*x+c)-3* 
b*c*(d*x+c)^2+b*(d*x+c)^3)+1/9*d^3/a/b^2*sum((3*_R1*b*c^2+a*d^3-b*c^3)/...
 
3.2.2.5 Fricas [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 0.30 (sec) , antiderivative size = 670, normalized size of antiderivative = 0.94 \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^3\right )^2} \, dx=-\frac {12 \, a d x \sin \left (d x + c\right ) + {\left ({\left (b x^{3} - \sqrt {3} {\left (i \, b x^{3} + i \, a\right )} + a\right )} \left (\frac {i \, a d^{3}}{b}\right )^{\frac {2}{3}} - {\left (b x^{3} + \sqrt {3} {\left (i \, b x^{3} + i \, a\right )} + a\right )} \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right )} {\rm Ei}\left (-i \, d x + \frac {1}{2} \, \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}} {\left (-i \, \sqrt {3} - 1\right )}\right ) e^{\left (\frac {1}{2} \, \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}} {\left (i \, \sqrt {3} + 1\right )} - i \, c\right )} + {\left ({\left (b x^{3} - \sqrt {3} {\left (i \, b x^{3} + i \, a\right )} + a\right )} \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {2}{3}} - {\left (b x^{3} + \sqrt {3} {\left (i \, b x^{3} + i \, a\right )} + a\right )} \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right )} {\rm Ei}\left (i \, d x + \frac {1}{2} \, \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}} {\left (-i \, \sqrt {3} - 1\right )}\right ) e^{\left (\frac {1}{2} \, \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}} {\left (i \, \sqrt {3} + 1\right )} + i \, c\right )} + {\left ({\left (b x^{3} - \sqrt {3} {\left (-i \, b x^{3} - i \, a\right )} + a\right )} \left (\frac {i \, a d^{3}}{b}\right )^{\frac {2}{3}} - {\left (b x^{3} + \sqrt {3} {\left (-i \, b x^{3} - i \, a\right )} + a\right )} \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right )} {\rm Ei}\left (-i \, d x + \frac {1}{2} \, \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}} {\left (i \, \sqrt {3} - 1\right )}\right ) e^{\left (\frac {1}{2} \, \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}} {\left (-i \, \sqrt {3} + 1\right )} - i \, c\right )} + {\left ({\left (b x^{3} - \sqrt {3} {\left (-i \, b x^{3} - i \, a\right )} + a\right )} \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {2}{3}} - {\left (b x^{3} + \sqrt {3} {\left (-i \, b x^{3} - i \, a\right )} + a\right )} \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right )} {\rm Ei}\left (i \, d x + \frac {1}{2} \, \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}} {\left (i \, \sqrt {3} - 1\right )}\right ) e^{\left (\frac {1}{2} \, \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}} {\left (-i \, \sqrt {3} + 1\right )} + i \, c\right )} - 2 \, {\left ({\left (b x^{3} + a\right )} \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {2}{3}} - {\left (b x^{3} + a\right )} \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right )} {\rm Ei}\left (i \, d x + \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right ) e^{\left (i \, c - \left (-\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right )} - 2 \, {\left ({\left (b x^{3} + a\right )} \left (\frac {i \, a d^{3}}{b}\right )^{\frac {2}{3}} - {\left (b x^{3} + a\right )} \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right )} {\rm Ei}\left (-i \, d x + \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right ) e^{\left (-i \, c - \left (\frac {i \, a d^{3}}{b}\right )^{\frac {1}{3}}\right )}}{36 \, {\left (a b^{2} d x^{3} + a^{2} b d\right )}} \]

input
integrate(x^3*sin(d*x+c)/(b*x^3+a)^2,x, algorithm="fricas")
 
output
-1/36*(12*a*d*x*sin(d*x + c) + ((b*x^3 - sqrt(3)*(I*b*x^3 + I*a) + a)*(I*a 
*d^3/b)^(2/3) - (b*x^3 + sqrt(3)*(I*b*x^3 + I*a) + a)*(I*a*d^3/b)^(1/3))*E 
i(-I*d*x + 1/2*(I*a*d^3/b)^(1/3)*(-I*sqrt(3) - 1))*e^(1/2*(I*a*d^3/b)^(1/3 
)*(I*sqrt(3) + 1) - I*c) + ((b*x^3 - sqrt(3)*(I*b*x^3 + I*a) + a)*(-I*a*d^ 
3/b)^(2/3) - (b*x^3 + sqrt(3)*(I*b*x^3 + I*a) + a)*(-I*a*d^3/b)^(1/3))*Ei( 
I*d*x + 1/2*(-I*a*d^3/b)^(1/3)*(-I*sqrt(3) - 1))*e^(1/2*(-I*a*d^3/b)^(1/3) 
*(I*sqrt(3) + 1) + I*c) + ((b*x^3 - sqrt(3)*(-I*b*x^3 - I*a) + a)*(I*a*d^3 
/b)^(2/3) - (b*x^3 + sqrt(3)*(-I*b*x^3 - I*a) + a)*(I*a*d^3/b)^(1/3))*Ei(- 
I*d*x + 1/2*(I*a*d^3/b)^(1/3)*(I*sqrt(3) - 1))*e^(1/2*(I*a*d^3/b)^(1/3)*(- 
I*sqrt(3) + 1) - I*c) + ((b*x^3 - sqrt(3)*(-I*b*x^3 - I*a) + a)*(-I*a*d^3/ 
b)^(2/3) - (b*x^3 + sqrt(3)*(-I*b*x^3 - I*a) + a)*(-I*a*d^3/b)^(1/3))*Ei(I 
*d*x + 1/2*(-I*a*d^3/b)^(1/3)*(I*sqrt(3) - 1))*e^(1/2*(-I*a*d^3/b)^(1/3)*( 
-I*sqrt(3) + 1) + I*c) - 2*((b*x^3 + a)*(-I*a*d^3/b)^(2/3) - (b*x^3 + a)*( 
-I*a*d^3/b)^(1/3))*Ei(I*d*x + (-I*a*d^3/b)^(1/3))*e^(I*c - (-I*a*d^3/b)^(1 
/3)) - 2*((b*x^3 + a)*(I*a*d^3/b)^(2/3) - (b*x^3 + a)*(I*a*d^3/b)^(1/3))*E 
i(-I*d*x + (I*a*d^3/b)^(1/3))*e^(-I*c - (I*a*d^3/b)^(1/3)))/(a*b^2*d*x^3 + 
 a^2*b*d)
 
3.2.2.6 Sympy [F(-1)]

Timed out. \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^3\right )^2} \, dx=\text {Timed out} \]

input
integrate(x**3*sin(d*x+c)/(b*x**3+a)**2,x)
 
output
Timed out
 
3.2.2.7 Maxima [F]

\[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^3\right )^2} \, dx=\int { \frac {x^{3} \sin \left (d x + c\right )}{{\left (b x^{3} + a\right )}^{2}} \,d x } \]

input
integrate(x^3*sin(d*x+c)/(b*x^3+a)^2,x, algorithm="maxima")
 
output
-1/2*(3*(cos(c)^2 + sin(c)^2)*d*x^2*sin(d*x + c) + ((d^2*x^3*cos(c) - 3*d* 
x^2*sin(c) - 12*x*cos(c))*cos(d*x + c)^2 + (d^2*x^3*cos(c) - 3*d*x^2*sin(c 
) - 12*x*cos(c))*sin(d*x + c)^2)*cos(d*x + 2*c) + ((cos(c)^2 + sin(c)^2)*d 
^2*x^3 - 12*(cos(c)^2 + sin(c)^2)*x)*cos(d*x + c) - 2*(((b^2*cos(c)^2 + b^ 
2*sin(c)^2)*d^3*x^6 + 2*(a*b*cos(c)^2 + a*b*sin(c)^2)*d^3*x^3 + (a^2*cos(c 
)^2 + a^2*sin(c)^2)*d^3)*cos(d*x + c)^2 + ((b^2*cos(c)^2 + b^2*sin(c)^2)*d 
^3*x^6 + 2*(a*b*cos(c)^2 + a*b*sin(c)^2)*d^3*x^3 + (a^2*cos(c)^2 + a^2*sin 
(c)^2)*d^3)*sin(d*x + c)^2)*integrate(3*(3*a*d*x*sin(d*x + c) + (a*d^2*x^2 
 + 10*b*x^3 - 2*a)*cos(d*x + c))/(b^3*d^3*x^9 + 3*a*b^2*d^3*x^6 + 3*a^2*b* 
d^3*x^3 + a^3*d^3), x) - 2*(((b^2*cos(c)^2 + b^2*sin(c)^2)*d^3*x^6 + 2*(a* 
b*cos(c)^2 + a*b*sin(c)^2)*d^3*x^3 + (a^2*cos(c)^2 + a^2*sin(c)^2)*d^3)*co 
s(d*x + c)^2 + ((b^2*cos(c)^2 + b^2*sin(c)^2)*d^3*x^6 + 2*(a*b*cos(c)^2 + 
a*b*sin(c)^2)*d^3*x^3 + (a^2*cos(c)^2 + a^2*sin(c)^2)*d^3)*sin(d*x + c)^2) 
*integrate(3*(3*a*d*x*sin(d*x + c) + (a*d^2*x^2 + 10*b*x^3 - 2*a)*cos(d*x 
+ c))/((b^3*d^3*x^9 + 3*a*b^2*d^3*x^6 + 3*a^2*b*d^3*x^3 + a^3*d^3)*cos(d*x 
 + c)^2 + (b^3*d^3*x^9 + 3*a*b^2*d^3*x^6 + 3*a^2*b*d^3*x^3 + a^3*d^3)*sin( 
d*x + c)^2), x) + ((d^2*x^3*sin(c) + 3*d*x^2*cos(c) - 12*x*sin(c))*cos(d*x 
 + c)^2 + (d^2*x^3*sin(c) + 3*d*x^2*cos(c) - 12*x*sin(c))*sin(d*x + c)^2)* 
sin(d*x + 2*c))/(((b^2*cos(c)^2 + b^2*sin(c)^2)*d^3*x^6 + 2*(a*b*cos(c)^2 
+ a*b*sin(c)^2)*d^3*x^3 + (a^2*cos(c)^2 + a^2*sin(c)^2)*d^3)*cos(d*x + ...
 
3.2.2.8 Giac [F]

\[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^3\right )^2} \, dx=\int { \frac {x^{3} \sin \left (d x + c\right )}{{\left (b x^{3} + a\right )}^{2}} \,d x } \]

input
integrate(x^3*sin(d*x+c)/(b*x^3+a)^2,x, algorithm="giac")
 
output
integrate(x^3*sin(d*x + c)/(b*x^3 + a)^2, x)
 
3.2.2.9 Mupad [F(-1)]

Timed out. \[ \int \frac {x^3 \sin (c+d x)}{\left (a+b x^3\right )^2} \, dx=\int \frac {x^3\,\sin \left (c+d\,x\right )}{{\left (b\,x^3+a\right )}^2} \,d x \]

input
int((x^3*sin(c + d*x))/(a + b*x^3)^2,x)
 
output
int((x^3*sin(c + d*x))/(a + b*x^3)^2, x)